Feshbach resonances in rubidium 87: Precision measurement and analysis 



A. 



(N 
O 
O 
(N 

o 

Q 



I 

o 
o 



(N 

o 

(N 
O 
-I— • 



X3 
O 

o 



X 



Marte\ T. Volz\ J. Schuster^ S. Durr\ G. Rempe\ E. G. M. van Kempen^, and B. J. Verhaar^ 
^ Max-Planck-Institut fiir Quantenoptik, Hans-Kopfermann-Str. 1, 85748 Garching, Germany 
^Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands 

(Dated: February 1, 2008) 

More than 40 Feshbach resonances in rubidium 87 are observed in the magnetic field range between 
0.5 and 1260 gauss for various spin mixtures in the lower hyperfine ground state. The Feshbach 
resonances are observed by monitoring the atom loss, and their positions are determined with an 
accuracy of 30 mG. In a detailed analysis, the resonances are identified and an improved set of model 
parameters for the rubidium interatomic potential is deduced. The elastic width of the broadest 
resonance at 1007 G is predicted to be significantly larger than the magnetic field resolution of the 
apparatus. This demonstrates the potential for applications based on tuning the scattering length. 

PACS numbers: 34.50.-s, 03.75.Fi, 32.80.Pj, 34.20.Cf 



A Feshbach resonance is an exciting tool for control- 
ling the atom-atom interaction in ultracold atomic gases. 
The elastic s-wave scattering length a can be tuned over 
orders of magnitude simply by applying a magnetic field. 
Feshbach resonances have been observed in various al- 
kali atoms H ||, |, |, |, 0]. They have been used to 
induce a controlled implosion of a Bose-Einstein conden- 
sate (BEC) to create a coherent superposition of an 
atomic BEC and a molecular state and to realize a 
bright soliton in a BEC |^ . Future applications could 
include experiments with the Mott-insulator phase tran- 
sition a Tonks gas effects beyond the mean- 
field theory and the creation of a molecular BEC 
p3[ . In addition, the binding energies of ro- vibrational 
molecular states close to the dissociation threshold can 
be determined with high accuracy from the position of 
Feshbach resonances leading to a precise knowledge of the 
interatomic potential [Q. A Feshbach resonance show- 
ing up in elastic collisions is often accompanied by strong 
changes in the inelastic collision properties |l6| . This 
offers a strategy for searching for new Feshbach reso- 
nances by monitoring the resulting atom loss. 

Surprisingly, no Feshbach resonance has been observed 
in *^Rb, which is the isotope used in most of today's 
BEC experiments. In 1995, a search with atoms in the 
\f,mf) — |1,— 1) state was carried out, but no reso- 
nances were found ||l7|| . Meanwhile various experiments 
[l9| , po| pT| greatly improved the knowledge of the 
Rb interatomic potential. Recent models based on this 
p2| , p3| are consistent with the observations of Ref. . 
But these models predict four Feshbach resonances in the 
1, 1) state. 

This letter reports the observation of more than 40 
Feshbach resonances in ^^Rb with most (but not all) 
atoms prepared in the |1,1) ground state. The theo- 
retical model is extended to various spin mixtures and 
to bound states with rotational quantum number I < 3. 
Thus, all except one of the resonances can be clearly iden- 
tified. Moreover, the measured position of one Feshbach 
resonance is used for an improved fit of the model pa- 



rameters. The relative deviation between the predicted 
and observed positions is 1.6 x 10"'^ (rms). The observed 
loss might be due to two- or three-body inelastic colli- 
sions, but should be purely three-body for those reso- 
nances which involve only the |1, 1) state. The broadest 
resonance at 1007 G offers the possibility to tune the 
scattering length and investigate its loss mechanism. 

The experiment is performed with a new set-up similar 
to our previous one , but with all relevant components 
significantly improved. In particular, atoms are captured 
in a vapor-cell magneto-optical trap (MOT) at a loading 
The atoms are transferred to a 



rate of 7 x 10^° s 



second MOT, in which 6 x 10^ atoms are accumulated by 
multiple transfer within 2 s. The atoms are then loaded 
into a loffe-Pritchard magnetic trap with a lifetime of 
170 s. The measured bias-field drift of less than 1 mC/h 
illustrates the excellent stability of the magnetic trap. 
After 26 s of evaporative cooling, a BEC with up to 3.6 x 
10^ atoms is formed in the |1, — 1) state. 

The |1, 1) state, in which the resonances are predicted, 
cannot be held in a magnetic trap. Therefore the atoms 
are now loaded into an optical dipole trap made of a 
single beam from an Yb:YAG laser at a wavelength of 
1030 nm. A laser power of 45 mW is focused to a waist of 
15 /im, resulting in measured trap frequencies of 930 Hz 
and 11 Hz, and an estimated trap depth of ks x 20 /xK. 

Once the atoms are in the optical trap, a radio- 
frequency field is used to transfer the atoms to the desired 

I, 1) state. With a Stern-Gerlach method, the fraction 
of atoms that end up in the |1,1) state is determined 
to be roughly 90%. Almost all other atoms are in the 

I I, 0) state. Next, a homogeneous magnetic field is ap- 
plied in order to observe a Feshbach resonance. The field 
is created using the compensation coils of the magnetic 
trap which are in near-perfect Helmholtz configuration. 
Up to 1760 A of current are stabilized with a home-built 
servo to a few ppm. The magnetic field is held at a fixed 
value for typically 50 ms and then quickly switched off. 
Although many Feshbach resonances are rapidly crossed 
when the magnetic field is turned on or off, no significant 
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TABLE I: Feshbach resonances in the |1, 1) dj) [1, 1) entrance 
channel. The experimentally observed positions Bexp are 
compared with the theoretical predictions Bth calculated at 
2 fiK. Also shown are the observed depths d, defined as the 
fraction of atoms lost during a 50-ms hold time. Note that 
the presence of atoms in spin states other than |1, 1) prevents 
d from reaching 100%. For comparison, theoretical results 
for the elastic widths A of the resonances are listed. The 
last columns list the quantum numbers of the (quasi) bound 
state that gives rise to the resonance. Some very weak I — 2 
resonances could not be detected experimentally. 



loss of atoms is observed from these rapid crossings. After 
switching off the magnetic field, the atoms are released 
from the optical trap, and 14 ms later an absorption im- 
age of the expanded cloud is taken. 

The search for Feshbach resonances was typically per- 
formed with a purely thermal cloud of 4 x 10^ atoms in 
the optical trap at a temperature of 2 /iK, corresponding 
to a peak density of 2 x 10^'' cni^^^. The magnetic-field 
range between 0.5 and 1260 G was scanned and 43 Fesh- 
bach resonances were found. The magnetic-field values 
of these resonances Bg^v are listed in Tabs. | and^. The 
fields were calibrated with 30 mG precision using mi- 
crowave spectroscopy in the vicinity of each resonance. 

In order to analyze these results theoretically, both 
continuum and bound-state calculations are carried out 
using an accurate description of the atomic interaction 
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between two Rb atoms 23 . An introduction to the the- 
ory of Feshbach resonances can be found in Refs. p6[ . 
The strong central part of the interaction, consisting of 



FIG. 1: / = Feshbach resonances in a coupled-channel 
calculation. Bound-state energies (solid lines) are shown as 
a function of magnetic field with quantum numbers (/i, f2)v' 
assigned at _B = 0. Additionally, dissociation threshold ener- 
gies (dotted lines) are shown for four different entrance chan- 
nels (mj",my^). A Feshbach resonance (•) occurs, when a 
bound state with quantum number mp crosses a dissociation 
threshold with the same mp. 



singlet and triplet potentials, conserves the orbital quan- 
tum numbers I, mi and the total-spin magnetic quantum 
number mp separately. (In bound-state spectroscopy, I 
is often called N.) In the ultracold regime, with this part 
of the interaction included, Feshbach resonances are ex- 
pected to occur only when an / = bound state with a 
certain value oi mp crosses the dissociation threshold of 
an entrance channel with the same mp. Figure |^ shows 
these crossings in the E-B plane as solid dots along vari- 
ous (my5_, w^^) thresholds. The resonance fields obtained 
from these bound-state calculations are listed in Tabs. || 
and H (/ = quantum number). The same field values 
also follow from the i?-dependence of the elastic S'-matrix 
element, as obtained in a continuum calculation. 

In Fig. |l|, the bound states are labelled (/i, f2)v' , with 
the vibrational quantum number v' — -1, -2, -3, ... 
counting from the corresponding (/i,/2) threshold. At 
-6 = 0, the exchange interaction couples the atomic spins 
fii /2 to a total molecular spin F and causes a split- 
ting between states with the same (/i, f2)v' but different 
F. In the presence of a strong external magnetic field, 
however, mfi and to/2 become good quantum numbers 
instead of F, while is always a good quantum num- 
ber. What constitutes a strong field in this sense depends 
on the size of the F-splitting at B = 0. For small \v'\, 
the F-splitting is hardly visible in Fig. 0. 

The much weaker spin-spin interaction Vss consists of 
the magnetic dipole-dipole interaction of the valence elec- 
trons together with a second-order spin-orbit term [ p7| . 
Due to its tensor form it breaks the spatial spherical sym- 
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TABLE II: Experimentally observed Feshbach resonances in 
other entrance channels \l,m^p^) ® |l,my^). 



metry and allows a redistribution between the angular 
momenta of the spin and spatial degrees of freedom so 
that the sum mi + mp is the only conserved quantum 
number. Vss admixes an Z = component in otherwise 
pure / — 2 bound states and therefore induces additional 
s-wave Feshbach resonances in the ultracold regime. Res- 
onances of this type have previously been observed in 
Cs 0, but up to now never in ^^Rb or *'^Rb. The re- 
sulting resonance fields are again listed in the tables {I 
= 2 resonances). For mixed species resonances, such as 
1, 0) 1 1, 1), the I = 1 partial wave can be populated in 
the entrance channel, thus opening up the possibility to 
observe resonances due to Z = 1 or Z = 3 bound states. 

The set of potential parameters used in this calculation 
is obtained as follows. The field value 911.74 G of one res- 
onance is added to the set of eight experimental data al- 
ready included in the analysis of Ref. . This particular 
resonance was chosen, because the corresponding bound 
state is a pure triplet state and the previous experimental 
constraints on the singlet potential are much stronger 
than those on the triplet potential. Then, a least-squares 
analysis is applied with the column A parameters of 
Tab. I in Ref. |Q as starting values. This leads to an ad- 
justed set of parameter values differing from the starting 
values by less than la; in atomic units: Cq = 4.707 x 10"^, 
Cg = 5.73 X 10^, Cio = 7.665 x 10^ (Ref. ||]), J = 
0.486 X 10-2, as(*^Rb) = +90.6, ari^'^Rh) = -1-98.96, 
?;Ds(^'^Rb) = 0.454, and vori^'^'Rh) = 0.4215. Coupled- 
channel calculations [^5[ |2^ based on these fine-tuned 
parameters then yield theoretical resonance fields Bth- 
The relative deviation of all these positions from the ob- 
servations is 1.6 X 10"'^ (rms). Compared to the four 
resonance positions predicted in |2^, this is an improve- 
ment of a factor of 6. Since only one resonance position 
was included in the fit, the excellent agreement with all 
other positions demonstrates the accuracy of the model. 



Note that the observed position of maximum loss might 
deviate from the pole of the elastic scattering length . 
However, such deviations are typically not larger than 
the elastic widths of the resonances; and they are small 
as discussed below. 

Some resonances are so close together that they cannot 
be identified merely from their positions. In these cases, a 
Stern-Gerlach method was used to experimentally deter- 
mine which m/ states incurred the strongest atom loss. 
Thus the entrance channel could be identified, in partic- 
ular when additionally varying the initial spin mixture. 
The entrance channel of the non-identified resonance at 
1236.73 G was also determined with this method. 

An interesting property of a Feshbach resonance - be- 
sides its position - is its strength. The strength of the 
elastic resonance is proportional to the field width A over 
which the scattering length has opposite sign as 
listed in Tab. |. The strength of the inelastic scatter- 
ing properties is quantified by the loss rate. The domi- 
nant loss for the resonances in Tab. | arises from inelastic 
three-body collisions. This is because for the parameters 
of the experiment, single-body loss is negligible; and since 
the |1, 1) entrance channel is the absolute ground state of 
atomic ^^Rb, inelastic two-body collisions cannot occur 
jlSf . (This is not the case for the resonances in Tab. ||.) 
The three-body loss is characterized by the coefficient 
in the rate equation N — —K^{'n?)N , where N is the 
atom number and n the density. The depth d listed in 
Tab. p is a non- linear, yet monotonic function of K^. In 
Tab. 1, one finds a clear trend that stronger resonances 
with larger A cause faster loss, i.e. larger d. This is plau- 
sible, although the theory of three-body losses is not yet 
fully understood, especially in the vicinity of Feshbach 
resonances |2^ . 

Interestingly, the widths A of the I — resonances 
in Tab. p| a re much smaller than for other alkali atoms 
iQ, H, |2l| . This is due to the approximate phase equal- 
ity of the waves reflected from the short-range singlet 
and triplet potentials, that is also responsible for other 
remarkable ^^Rb phenomena: the coexistence of conden- 
sates in different hyperfine states ||T^ and the smallness 
of the fountain-clock frequency shift p^ . 

All resonances display a nearly symmetric loss feature. 
A Gaussian is fit to the atom loss at a 50-ms hold time 
in order to extract the depth and the width of the atom 
loss. The Gaussian was chosen for convenience and be- 
cause it fits well to the data. The obtained rms-widths 
are between 20 mG and 100 mC for almost all resonances. 
For small d, this width is identical to the width of the 
resonance in i^a (if two-body loss is absent). For large 
d, however, the non- linear dependence of d on broad- 
ens the observed width at a given hold time, somewhat 
similar to saturation broadening of spectral lines. In ad- 
dition, the finite temperature of the cloud gives rise to 
a broadening of typically 20 mC and leads to a shift of 
similar size. The observed rms-width of 24(4) mC at 
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FIG. 2: Shape of the broadest Feshbach resonance. The 
number of atoms (o) remaining in the trap after a 50-ms hold 
time is displayed as a function of magnetic field. The three- 
body loss rate coefficient Ks is also shown (•). 



965.96 G sets an experimental upper limit to thermal 
and technical broadening. 

The broadest resonance is centered at 1007 G (see 
Fig. H). The theoretical prediction for its elastic width, 
A = 170(30) mG, is large compared to the upper bound 
on the experimental broadening of 24(4) mG (see above). 
This demonstrates the potential for a controlled variation 
of the scattering length with the present set-up. 

For various field values, the decay of the atom num- 
ber was also measured as a function of time. If one 
assumes that the loss can be described by a rate equa- 
tion, three-body loss dominates (see above). Values of 
K3 determined from a fit are shown in Fig. ^. Note 
that the value K3 = 3.2(1.6) x lO^^g cm^s^^ obtained 
away from the resonance is consistent with the measured 
value K3 = 4.3(1.8) x lO'^^ cm^s'^ for the |1, -1) state 
p0| . The absence of two-body loss makes this resonance 
an ideal candidate for testing theories of three-body loss 
pgf . The values of K3 in Fig. ^ are much smaller than 
those for ^^Rb |16[, where exciting experiments have been 
performed ||, | ] . Three-body loss will therefore not be a 
substantial problem for applications. 

All data obtained here fit well to a three-body decay, 
but the initial atom number extrapolated from the fit 
is a factor of up to 2 lower at the 1007 G resonance as 
compared to the off-resonance value. Since no data were 
taken for hold times shorter than 10 ms, this suggests the 
existence of an additional loss mechanism acting on much 
faster timescales. This could be molecule formation and 
dissociation during the first 10 ms of hold time Note, 
however, that here loss from a purely thermal cloud is 
observed, while loss from a BEG is discussed in Ref . Q . 

To summarize, more than 40 Feshbach resonances have 
been observed in ®^Rb. All except one were identified by 
theory. After including one of the observed resonances 
in fitting the model parameters, theory and experiment 
are in excellent agreement. The magnetic field control is 



so accurate that it should be possible to resolve changes 
in the elastic scattering length on the broadest resonance 
at 1007 G. 
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